Anomalous elasticity of nematic elastomers 
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We study the anomalous elasticity of nematic elastomers by employing the powers of renormalized 
field theory. Using general arguments of symmetry and relevance, we introduce a minimal Landau- 
Ginzburg- Wilson elastic energy for nematic elastomers. Performing a diagrammatic low temperature 
expansion, we analyze the fluctuations of the displacement fields at and below the upper critical 
dimension 3. Our analysis reveals an anomaly of certain elastic moduli in the sense that they depend 
on the length scale. In d = 3 this dependence is logarithmic and below d = 3 it is of power law 
type with anomalous scaling exponents. One of the 4 relevant shear moduli vanishes at long length 
scales whereas the only relevant bending modulus diverges. 

PACS numbers: 61.41. +e, 64.60.Fr, 64.60.Ak 



Nematic elastomers jl], ||, || are elastic media with the 
macroscopic symmetry properties of nematic liquid crys- 
tals [0. In addition to the degrees of freedom of ordi- 
nary elastomeric materials, liquid-crystalline elastomers 
possess an internal, orientational degree of freedom. For 
moderate crosslinking, the existence of the rubbery net- 
work has apparently little impact on the liquid crystalline 
phase behavior. However, because liquid-crystalline elas- 
tomers cannot flow, they have mechanical properties that 
differ significantly from standard liquid crystals. 

Nematic elastomers have unique properties that make 
them candidates for device applications. Temperature 
change JfJ or illumination Q can change the orientational 
order and cause the elastomer to extend or contract as 
much as 400% M. This qualifies nematic elastomers as 
contestants for use in artificial muscles ||, [I. Moreover, 
nematic elastomers display a soft elasticity flO], [Tl|| char- 
acterized by vanishing shear stresses for a range of lon- 
gitudinal strains applied perpendicular to the uniaxial 
direction. 

The elastic energy of nematic elastomers was first stud- 
ied by Golubovic and Lubensky (GL) |To| . Normally, 
the elastic energy of a uniaxial elastic medium features 
five independent elastic constants. However, when a 
nematic elastomer is synthesized by crosslinking in the 
isotropic phase, the uniaxial state arises from a spon- 
taneous symmetry breaking at the isotropic to nematic 
transition. This has the important consequence that the 
elastic constant associated with shears in the plane con- 
taining the anisotropy axis vanishes. The work of GL in- 
dicated that the upper critical dimension d c for nematic 
elastomers is 3 or 5 in the absence or presence of ran- 
dom stresses, respectively. In each case, one can expect 
anomalous elasticity at or below the respective d c due to 
Grinstein-Pelcovits |Q type renormalizations of the re- 
maining elastic constants. Here, we study the anomalous 
elasticity of nematic elastomers without random stresses 
by using the powerful methods of field theory [[15] aug- 



mented by a renormalization group (RG) analysis. Our 
aim is to derive the anomalous scale dependence of the 
elastic moduli in and below 3 dimensions. Parallel to 
our work, Xing and Radzihovsky independently studied 
nematic elastomers without jl3j and with [ pd| random 
stresses. Though obtained by an entirely different for- 
malism their final results p3| completely agree with ours. 

For our discussion, we first must devise an appropri- 
ate minimal model for nematic elastomers in form of 
a field theoretic Hamiltonian. Following GL, we use 
the Lagrangian formulation of elasticity ]l6| ]. In this 
formulation the mass points of the equilibrium undis- 
torted medium are labeled by their position vectors x 
in cZ-dimensional (reference) space. When the medium is 
distorted, a mass point originally at x is mapped to a 
new point R(x) in d-dimcnsional (target) space. Since 
R(x) = x when there is no distortion, it is customary 
to introduce the phonon variable u(x) = R(x) — x that 
measures the deviation of R(x) from x. The energy of 
the distorted state relative to the reference state depends 
only on the relative amount of stretching <iR 2 — dx 2 = 
2uijdxidxj, where the Einstein convention is understood 
and where the uy = ^ (diUj + djUi + diUkdjUk) with 
i,j,k = l,...,d are the components of the familiar non- 
linear Lagrangian strain tensor u. By construction u is 
invariant under arbitrary rotations in target space. This 
feature makes the Lagrangian strain tensor an adequate 
variable for formulating elastic energies because all elastic 
media are rotationally invariant in target space. This in- 
variance is easy to understand: in the absence of aligning 
fields different physical orientations of the same sample 
have the same energy. 

We are interested in nematic elastomers crosslinked in 
the isotropic phase whose uniaxial symmetry arises via 
a spontaneous breaking of the isotropic symmetry. Since 
the spontaneous symmetry axis can point in any direc- 
tion, our elastic energy has to be rotationally invariant 
not only in the target but also in the reference space. 
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Both invariances are taken into account by writing the 
stretching energy as 

H = J d d x^(trgf + fitrg 2 + A^trg) 3 

+ A 2 tr^tr^ 2 + A 3 tr^ 3 + Bi(tr^) 4 + B 2 (trg 1 ) 2 
+ B 3 (tr^) 2 tr^ 2 + B 4 trjitrj| 3 + B 5 tr^ 4 | . (1) 

The first two expansion coefficients are known as the 
Lame coefficients. Of course, terms of higher than fourth 
order are allowed by the symmetries of the system. How- 
ever, these higher order terms turn out to be irrelevant 
in the RG sense and are hence neglected. 

Suppose that the spontaneously uniaxially ordered 
elastomer is described in equilibrium by an equilibrium 
strain tensor u^. Without loss of generality we may as- 
sume that the anisotropy axis lies in the = (0, . . . , 1) 
direction and that u Q is a diagonal matrix with the di- 
agonal elements uoaa — uqj_ and uqm = Uq\\- Here and 
in the following we use the convention that indices from 
the beginning of the alphabet (a, b, c) run from f to 
d\_ = d — 1. To describe deviations from the equilibrium 
configuration, we introduce the relative strain w = u — u Q 
and to expand Ji in terms of w. By dropping terms that 
depend only on u Q we find 

^ = / d<ix { aiWdd + a2Waa + b l W dd + hWddWaa 

+ b 3W 2 aa + b A w 2 ab + b 5 w 2 ad + dw dd wl d + c 2 w aa w 2 bd 
+ c 3 w ab w ad Wbd + diwl d w 2 d } , (2) 

where J d d x — J d d± x J dx d with dj_ = d — 1. In (Q) we 
have discarded terms that turn out to be irrelevant. The 
new coefficients ai, a 2l bi, and so on, depend on the old 
coefficients A, fj,, and so forth, as well as on u || and uq±. 

Exploiting the rotational invariance in reference space 
we derived a set Ward identities relating the vertex func- 
tions implicit in (^). At zero-loop or mean- field level 
these Ward identities correspond to relations between the 
elastic constants in (^), 

ai — a 2 — s b 5 — , (3a) 
b 2 - 2b 3 - s c 2 = , b 2 + b 5 - 26i + s c x = ,(3b) 
65 — 264 — s C3 = , ci — C2 — C3 — 2s di = ,(3c) 

where s is an abbreviation for u ii — uq±. Since w de- 
scribes deviations from the equilibrium m q its thermal av- 
erage (w) has to vanish. At zero temperature, where the 
mean-field approximation becomes exact, this means that 
the coefficients of the linear terms in (Q) must be zero. 
Equation ( |3a| ) then leads to the observation that 65 = 
for u || 7^ uq±. At finite temperatures thermal fluctua- 
tions become important and loop corrections renormalizc 
the elastic constants including oi, a 2 and 65. For (w) to 
vanish the renormalized versions of a\ and a 2 have to 



satisfy equations of state to which a\ = and a 2 = 
are the mean- field approximations. In the following we 
will assume that we have chosen a± and a 2 appropriately 
so that their respective equations of state are satisfied. 
In other words: we assume that we expand about the 
true equilibrium state. Then, the Ward identity general- 
izing (|3a|) guarantees that the renormalized 65 vanishes 
for u || =/= uqj_. The vanishing of the elastic constant 65 
has the striking consequence that nematic elastomers are 
soft with respect to certain deformations. This softness 
can be easily understood by looking at the terms of (|^) 
that are leading for small deformations. Rewriting these 
terms in Fourier space, one sees that that the stretching 
energy cost is zero for phonon displacements u(q) per- 
pendicular to with momentum q parallel to and for 
u(q) parallel to with q perpendicular to e^. 

For vanishing 01, a 2 and 65 one can rewrite 7i after a 
suitable rescaling of x d and u d as 

H = J d d X^^-V 2 dd + y (y\u d ) 2 + C 2 V dd V aa 

+ ^ 2 aa + C iV 2 ab }, (4) 

with Vab = \{d a Ub + d b u a - d a u d d b u d ) and v dd = d d u d + 
\d a u d dbU d . Here, we exclusively retained terms that rel- 
evant in the RG sense. The new elastic constants C... 
are proportional to the 6... in (S). Note that we have in- 
corporated a relevant bending term featuring a bending 
modulus K . Bending is important in nematic elastomers 
since the pure stretching energy vanishes for soft defor- 
mations. Note also that Tt resembles the form of the 
stretching energy originally discussed by GL. However, 
we have to keep in mind that the strain v is different 
from u. Ti. reduces to the elastic energy by GL if one 
neglects terms of higher than second order in u. 

In principle we could use H as it stands in Eq. (Q) for 
our RG analysis. We find it convenient, however, to re- 
duce the number of constants featured in the statistical 
weight exp(— Ti/T) at the onset of our calculation. To 
this end, we rescale T \J K 3 j C4T , x d — * \fC4~l Kx d , 
u d — > \J KjC^u d and u a — * {KjC^)u a . After the rescal- 
ing, there remain 3 (dimensionless) couplings u> = C1/C4, 
g = C2/C4, and / = C3/C4 in our Hamiltonian. Note 
that we choose to keep T explicit and to use it as the 
expansion parameter for our field theoretic perturbation 
expansion. 

For our RG analysis, two scaling symmetries are of 
great importance. First, under a global rescaling of the 
coordinates x a —> pT 1 x a and x d — > \x~ 2 x d with hav- 
ing units of length, we find a scaling invariant theory 
provided that u a — > (J^u a and T — > fi 3 ~ d T. Thus, the 
naive dimension of u a is I and that of T is e = 3 — d. 
The field u d and the remaining parameters in 7i have 
vanishing naive dimensions. Above d — 3 dimensions, 
the naive dimension of T is negative, and T is irrele- 
vant whereas it is relevant below d — 3. Hence, we 
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identify d c — 3 as the upper critical dimension of ne- 
matic elastomers. Second, TL is invariant under the trans- 
formation u a {x ci x d ) -> u a (x c - 6 c x dl x d ) + 9 a u d (x ci x d ) 
and u d {x c , x d ) -> u d (ir c - 6 c x d , x d ) + 8 a x a provided that 
the #'s are small. Note that this transformation mixes 
the longitudinal and the transversal fields (mixing invari- 
ance). It can be viewed as a rudiment of the rotational 
invariance of the original theory in target space. This 
mixing invariance is valuable because it guarantees Ward 
identities that reduce the number of vertex functions to 
be calculated in perturbation theory. 

Now we study fluctuations of the fields u a (pc) and u d (x) 
by using perturbation theory augmented by renormaliza- 
tion group methods. As usual, we analyze vertex func- 
tions that require renormalization due to the presence 
of ultraviolet (UV) divergences in Feynman diagrams. 
To avoid infrared (IR) singularities in the Feynman dia- 
grams, we supplement our Hamiltonian with a mass term 
J J d d xu 2 d . At the appropriate stage of the calculations 
we then sent r to zero to recover the original situation. 

In order to set up a diagrammatic perturbation ex- 
pansion we have to determine its constituting elements. 
First, we have a Gaussian propagator that has the form 
of a d x d matrix. Its matrix elements are 
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Once the equations of state are satisfied and the 
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(5b) 
(5c) 



where we have used the shorthand notations A — 
w(J + 2) - g 2 , B = f + 2, A = u(f + 1) - g 2 , 
and B = / + 1. Second, our diagrammatic expan- 
sion features the 4 vertices i(u> — g)/(2T) <ld% % \ 
i(g - f)/(2T) qa 1] q ( b 2) qf\ -ih/T $\ and -(co - 

2g + f + 2)/(8T) q a 1] qi 2) qf ] q { b ] ■ It is understood that the 
sum of the momenta at each vertex has to be zero. 

Next we need to determine which of the vertex func- 
tions r^'"' with N external u a -legs and M ex- 
ternal Mrf-legs are superficially UV divergent. Analyzing 
their topology, we find that the superficial degree of di- 
vergence (5 of our diagrams is given at the upper critical 
dimension by 8 — 4 — N — 2Du — D±, where D» (D±) 
is the number of longitudinal (transversal) derivatives on 
the external legs. Thus, the vertex functions containing 



(0,1) 



r (i,o) 

1 a , 



,(0,2) 
dd 



superficially divergent diagrams are T d 
p(l,l) p(2,0) p(0,3) p(l,2) , r (o,4) A|] , 
1 ad I 1 ab > 1 ddd > 1 add > &nCX 1 dddd' A11 tneSe v ^ n ^ x 

functions have to be taken into account in the renormal- 
ization procedure. By virtue of the mixing invariance, 
however, there exist several relations between the ver- 
tex functions in form of Ward identities. Due to these 
Ward identities it is sufficient for our purposes to ac- 



2-point functions are renormalized, the Ward identities 
guarantee that the remaining vertex functions are cured 
of their UV divergences. 

We calculate the two-point vertex functions to one-loop 
order using dimensional regularization. In dimensional 
regularization the UV divergences manifest themself in 
the form of poles in the deviation e = 3 — d from the 
upper critical dimension. We eliminate the e poles from 
the vertex functions by employing the renormalization 
scheme 



ry-l/2 

Xd -> Xd = Z 1 x d , 
u a ^ u a = Zu a , 
lu — ► lu = Z 1 Z^p 1 Z^jUJ , g 

f - / = z- 1 z T 1 z f f , 



u d — > u d = Z 1/2 u d , (6a) 
T^f= Z 1/2 Z^ 1 ^ti&o) 
g = Z- 1 Z T 1 Z g g, (6c) 
(6d) 



where the ° indicates unrenormalized quantities. The sim- 
plest way of determining the renormalization Z-factors is 
minimal subtraction. In this procedure the Z-factors are 
chosen so that they solely cancel the e poles and other- 
wise leave the vertex functions unchanged. Our Z-factors 
are of the structure Z, .. (t, u, g, /) = 1 + tX,..(w,g, f)/e, 
where the X._(oj,g, f) are rational functions. 

Next, we set up a Gell-Mann-Low RG equation (RGE) 
for the vertex functions by exploiting the fact that the 
unrenormalized theory has to be independent of the ar- 
bitrary length scale introduced by renormalization. 
At this stage we prefer to switch from the parameters w, 
g and / to K = g/u> = C2/C1, p = f/uj = C3/C1 and 
(7 = 1/oj = C4/C1. This step turns out to be helpful in 
studying the RG flow. Our RGE reads 

{2? M -(N + M/2) 7} T (Ar ' M) ({q ± , q d } ; t, «, p, a, M ) = , 

where 2? M = pd^ - \q d d qd + (3 t d t + (3 K d K + (3 p d p + (3 a d a , 
ft... = M^mIo-- an d 7— = A*^mIo m Z. .. 0- To solve the 
RGE we employ the method of characteristics. We intro- 
duce a flow parameter £ and look for functions p(£), Z(t), 
i(£), k(£), p(£) and a(£) determined by a set of charac- 
teristic equations. These equations specify how the pa- 
rameters transform if we change the momentum scale p 
according to p, — * p(£) — tp. In particular, we have 



£d e t 

Idt R 
id t p 

Idea 



-te + P[-3R 2 -12ka 
p{i + 37a)]/[CTfiA] , 
i(R - l)(p + a- R 2 )4n/A, 
t[2R 2 (l -p) + 2f - 4Ra + a 2 
p(4a--2)]2fi/A, 
ia(-2 + 4k - 2R 2 + a) 2fl/A , 



2ct(3 + 19a) 



(7a) 
(7b) 

(7c) 
(7d) 



where O = + p/a and A = 64 / K^/~p~- 
First, we consider the case d < 3. 



tually calculate the 2-point functions T 



(0,2) p(l,l) 



dd 



and 



2(7 — R . 
< 3. Being interested 
in the IR behavior of the theory, we study the limit 
I — > 0. We find that the quadruple of coupling con- 
stants (i(£), R(£),p(£), <j(£)) flows, in a manner such that 
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k(£) — k* ~ p(£) - p* ~ cr(£) for small to an IR sta- 
ble fixed point (t*,K*,p*,a*) = (flv^Tre, 1,1,0). Note 
that this fixed point implies the existence of the following 
universal ratios: C2/C1 = 1, C3/C1 = f and C4/C1 = 0. 
These should be conveniently measurable in experiments. 
In addition to the stable fixed point, there is the unstable 
Gaussian fixed point t* = and there are two unstable 
fixed lines. Solving the RGE near the stable fixed point, 
in conjunction with dimensional analysis to account for 
the naive dimensions, reveals the scaling behavior of the 
vertex functions, 

T (N,M) ({ q ^ qd } ;t ,K,p,a,p) 

_ £-(d+l)+dN+(d+l)M Z (£y(N+l)/2 



In summary, we have derived a field theoretic model 
for nematic elastomers and explored, on this basis, their 
fascinating anomalous elasticity in and below 3 dimen- 
sions. Our work underscores the intriguing properties 
of nematic elastomers and calls for experimental verifi- 
cation. We expect that our general approach is fruitful 
also for other liquid crystalline elastomers and that cor- 
responding studies will appear in the future. 
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00-95631(TCL). We are grateful to H. K. Janssen, J. 
Toner, X. Xing, and L. Radzihovsky for fruitful discus- 
sion and to the latter two for pointing out an inconsis- 
tency in a preliminary version of this paper. 



where Z{£) ~ £~™ £ and a{£) ~ m £ . 

Now, the scaling behavior of the elastic constants can 
be easily deduced from the two point functions. Choosing 
£ = |q_|_| and switching back to the original variables 
featured in the elastic energy (^) we find for example 
that 

ri°/ } (qx, to) = T- 1 {Ch q 2 d + K(q ± ) |qj 4 } (9) 

with a normal (constant) shear modulus C\ and a mo- 
mentum depending bending modulus if(q^) ~ |q_i_| ~' ,A 
governed by the anomaly exponent r\K = §§£■ Note that 
K diverges at large length scales. The scaling behavior of 
the remaining elastic constants can be derived by analo- 
gous means. We find that C2 and C3 are normal whereas 
C,i(q_i_) ~ IqjJ 1 * 4 is anomalous with 774 = j4e, i.e., C4 
vanishes a large length scales. We mention that (|J) docs 
not correctly describe T^' - 1 at q_|_ = because we omit- 
ted a dangerous irrelevant bending modulus. This intri- 
cacy, however, does not effect our results for the relevant 
clastic moduli and will be discussed elsewhere. 

Since e vanishes in d = 3 the solutions to the character- 
istic equations are no longer of power law type. The flow 
of the temperature, for example, is described at leading 
order by t(£)/t = [1 - In^)]" 1 where t = t(l). Simi- 
larly, we have Z{£) ~ [t{£)/t]~^ and a(£) ~ [i(£)/t}& at 
leading order. Equation m) now implies that 



^(q±) 

C 4 (q ± ) 



7V6t 
64tt 

64tt 



m(/V|qJ) 



ln(/V|qJ) 



38/59 



-4/59 



(10a) 



(10b) 



Of course, C\, C2, and C3 are normal as they are in d < 3 
and K and C4 diverge and vanish, respectively, at long 
length scales. Note that Eqs. ( J10| ) imply the existence of 
a nonlinear length scale £ t = // exp( 7 ^ f ). 
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